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SUMMARY 


The linearized drag integral for bodies of revolution at supersonic 
speeds is presented in a double-integral form which is not based on 
slender-body approximations but which reduces to the usual slender-body 
expression in the proper limit. With the aid of a suitably chosen aux- 
iliary condition, the minimum-external-vave-drag problem is solved for 
a transition section connecting two semi-infinite cylinders. The pro- 
jectile tip is a special case and is compared with the Von K£rman pro- 
jectile tip. Calculations are presented which indicate that the method 
of analysis gives good first-order results in the moderate supersonic 
range. 


INTRODUCTION 


In maki ng the slender-body approximation to the linearized 
supersonic- flow theory for bodies of revolution, a basic approximation 
leads to replacing the axial source distribution with the cross-sectional- 
area derivative. Slender-body theory therefore becomes linear in the 
superposition sense for cross-sectional areas as well as for sources 
or fields in contrast to linear supersonic-flow theory which is linear 
in the superposition sense for sources or fields but not for areas. 

Making the slender-body approximation, however, eliminates a large part 
of the Mach number dependence of the results. Lighthill (ref. l) has 
shown that this basic approximation, for sufficiently smooth bodies, 
has the same mathematical order of accuracy as the linearized supersonic- 
flow equation. Ward (ref. 2) has extended the generality of slender- 
body theory by presenting a drag expression which is valid for a body 
with a finite slope at the base. Tight hi 11 (ref. 3) has, at the price 
of a large increase in complexity, modified slender-body theory to 
include area-derivative discontinuities at a finite number of points. 

Tn 1935 Von Karman (ref. 4) determined the minimum-wave-drag pro- 
jectile tip. Later Sears (ref. 5) and Haack (ref. 6} determined minimum- 
wave-drag shapes for projectile tips and closed bodies of revolution 
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subject to various combinations of auxiliary conditions of constant 
length, constant caliber, and constant volume. The Von Karman tip and 
the Sears-Haack bodies are based on slender-body theory. In contrast 
to the minimum-drag bodies of given volume, which are accepted as 
reasonably good first-order results, the minimum-drag projectile tips 
and open-nosed bodies which have been derived from slender-body theory 
are subject to more severe restrictions on the range and conditions for 
validity. Busemann (ref. 7 ) pointed out this fact for the Von Karman 
projectile tip in 1941. In the late 1930 's Ferrari made studies of the 
mini mum-drag projectile tip and length- caliber body problems. Recently, 
Ferrari (refs. 8 and 9) bas considered the length-caliber body and 
ducted-body minimum-drag problems on the basis of linear theory without 
resorting to slender-body approximations. He obtains source-distribution 
functions which involve elliptic integrals that have to be treated 
numerically. 

This paper presents the llnearized-drag integral for bodies of 
revolution in a double-integral form which is not based on slender- 
body approximations but which reduces to Ward's drag expression in the 
proper limit. With the aid of a suitably chosen auxiliary condition, 
the minimum-extemal-wave-drag problem for a transition section con- 
necting two semi- infinite cylinders is solved. The source distribution, 
the minimum-wave drag, the slopes at the ends of the section, and the 
radius at one intermediate point are obtained in terms of elementary 
functions. The entire shape Is given in an integral form amenable to 
numerical evaluation. 
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SYMBOLS 

constant factor in source-distribution function, ft /sec 
drag coefficient 
drag, lb 

source-distribution function 

derivative of f with respect to its argument 

free-stream Ifech number 

radius in cylindrical coordinates 

general radius of point on body 

radius of downstream cylinder 



NACA TN 3189 


3 


r 2 

A'(x),A"(x) 

u 

v r 

v x 

X 

Xq 

X 


radius of upstream cylinder 

first and second derivatives of cross-sectional area 
free-stream velocity, ft/sec 
radial velocity, ft/sec 
axial velocity, ft /sec 

axial distance in cylindrical coordinates 
position of first disturbance on control surface 
axial distance for general point on body 


0 = i/m 2 - 1 

Sf ' 

P 

A 

£0^ ^o 
0 


variation in derivative of source-distribution function 

free-stream density, slugs/cu ft 

arbitrary constant 

dummy variables of integration 

position on axis where sources begin 

disturbance potential 


All distances are made dimensionless by measuring them in the units 
of length of the transition section. 


ANALYSIS 
Drag Integral 


The linearized-flow equation for the case of axial symmetry is 
(ref. 4) 


£20 + 1 80 

£>r 2 r Sr 



( 1 ) 
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where $ is the disturbance potential, p^ = - 1, and M is the 

free-stream Mach number. The general solution of equation (l) is 
(ref. 10) 


0(x,r) = f f(x - j3r cosh u)du = - f P — f ^ d ! 


/(x - §) 2 - p 2 r 2 


( 2 ) 


The corresponding disturbance velocities are 


I* - v x (x,r) - r° f’(x - Sr cosh u)'4u = . f'tiHs (J) 

831 J " J -“ /(x - f) 2 - 


and 


^ = v r (x,r) a .f°f(x. p r cosh u)p cosh udu = i f' < U ■“ 

Sr /(x -if - p 2 r 2 


( 4 ) 


The general solution is interpreted physically as the disturbance 
potential of a distribution of axial sources. The source-distribution 
function f(|) is fixed by satisfying the boundary condition that the 
body surface be a stream surface. 

Applying the momentum and continuity theorems to the fluid within 
a cylindrical surface of radius equal to the radius of the downstream 
cylinder (hereafter referred to as the control surface) and using 

the fact that free-stream conditions prevail at the two ends of the 
cylinder, x = » and x = -«>, yield the usual expression for the drag 
(ref. k) 


D = -2rt 



PRlv x (x,Ri)v r (x,R 1 )dx 


( 5 ) 


If a position x = L on the control surface is chosen so that, for 
x ^ L, either v r or v x vanishes on the control surface, and if 
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Xq is the position of the first disturbance on the control surface, 
equation ( 5 ) reduces to 


D = -2rt J' pR 1 v x (x,R 1 )v r (x,R 1 )dx 


( 6 ) 


Substituting expressions ( 3 ) and (4) for the velocities into equation (6), 
■where £0 = T k> = *0 " the position on the axis where the sources 

begin, yields 


D = 2jtp 


r 

nx-pRi f'( 5 )(x _ !) d g 

r x ~ m 1 

x o 

Jo j/(x - |) 2 - 

Jo ^x - T ]) 2 - p2R L 2 


dx 


or 


D = 2itp 


f 


'L 

Xo 


rx-gRp rx-pRi f '(£)f ' (n)(x - |)d£ dr] 

£ o Vo j/(x - |) 2 - p% L 2 /(x - ii) 2 - p%i 2 


(7) 

dx (8) 


The basic slender-body approximation reduces the integral over | 
in equation ( 7 ) to f(x). The use of equation (4) for small values of r 

and the boundary condition ^ = ^jr- gives f(x) proportional to the 

axial derivative of the cross-sectional area. The slender-body approxi- 
mation thus reduces the integral of equation (8) to a double integral 
(over x and T}). In principle, the integral in equation (8) may be 
reduced to a double integral over £ and rj by interchanging the order 
of integration and performing the integration over x. When this process 
is attempted in equation (8), however, the integration over x involves 
elliptic integrals and no essential advantage is gained. 

However, the drag integral may be reduced to a double integral in 
the following manner. Equation (8) is rewritten with the dumny vari- 
ables | and T] interchanged. Addition of the interchanged and 
original expressions and division by 2 yield 


D 



f ’(|)f ' (ti)(x - | + x - Ti)d| di) 

/(x - i) 2 - (3%i 2 |/(x - T)) 2 - 0%i 2 


(9) 



Interchanging the order of integration yields 
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*JLfl 
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v* 

11 

P 


Equation (12) gives the drag of an axial distribution of sources f(g) which begins at 

and is subject to the condition that there is no momentum flow through the control cylinder of 
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In the slender-body approximation, equation (12) reduces to Ward's 
drag expression as follows. Use of the expression for the radial 
velocity (eq_. (if) ) for small values of r and the boundary condition 

dr v r , 

~ = ~ give 
dx U 


f(x) = ~ A'(x) 


( 13 ) 


where A'(x) 
small values 
integrals in 

by 


is the axial derivative of the cross-sectional area. For 
of R]_, pRj_ may he neglected in the limits of the 
equation (12) and the cosh“l factor may he approximated 


cosh 


-1 


(L - £)(L - rQ - p^i 2 

pRl(ft - n) 


= log [7l - £)(L - nfj - 
log ji - T I - log 


Cl 1 !-) 


Simple manipulation, based on the assumption that A’(|) 1 b continuous 
and A'(| 0 ) = 0, gives the Ward drag formula 


D 


“ J J a" ( §)a" (7])log J | - T] jdg dT] + 

to Tlo 

~ A'(L) f A"(|)log(L - |}d| - HH!jA’(Lfj 2 log ^ 

2rt Jt kft 2 


( 15 ) 


Minimum-Drag Problem 

The problem of determining the shape of a transition section con- 
necting two semi- inf ini te cylinders for minimum external wave drag may 
be formulated in the following way. (See fig. 1.) The upstream and 
downstream semi-infinite cylinders are of radius Rg and R^, respec- 
tively. The length of the transition section Is taken as the unit in 
which all distances are measured. The cylinder of radius Rq is chosen 
as the control surface. 
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Figure 1 is drawn for the case R 2 < Rq. However, the statement 
Of the problem is the same as for Rg > Ri since only the appropriate 
changes in sign occur. Equation (12) for the drag becomes 


D = up 


AL-3Rq rl-pRq 

/ / f'U)f (ii}cosh“ 1 

(1 - 

|)(1 - T|) - pSRq 2 

; -pR 2 J -3R 2 


PRqU - ri) 


d| dq 


(16) 

The auxiliary condition for the minimum-wave-drag problem is the 
requirement that the points x = 1, r = Rq and x = 0, r = R 2 be 
on the same stream surface. In terms of mass flow this condition 
requires that the free-stream mass flow between the control cylinder 
of radius Rq and the upstream cylinder of radius R 2 equals the maBB 
flow out of the control cylinder between x = p(Rq - R 2 ) and x = 1; 
that is, 


pUlt(Rq 2 - R 2 2 ) 


/ 

J p(Rq-R 2 ) 


2rtpRqv r (x,Rq)dx = Constant 


or 


(17) 


f 

p (Rq-R 2 ) 



t'UKx -itel . ^ = U( Rl 2 _ R2 2) . ConBtant 

/(* . I) 2 . p2ftq2 


(18) 


Equation (18) is a particularly convenient form of the auxiliary 
condition in that the source-distribution function is the only unknown 
involved. A direct integration of the usual boundary condition would 
lead to an auxiliary condition involving, in addition, the unknown 
shape of the section. Without further approximations, the minimum- 
drag problem would be unmanageable. It appears impossible to obtain 
an auxiliary condition of constant volume involving the source-distribu- 
tion function as the only unknown. 

Interchanging the order of integration in the left member of equa- 
tion (l8) and integrating over x yield 


f f '(I) - |} 2 - p^q^i = |(Ri 2 - R 2 2 ) = Constant 

w _ o*D _ 


( 19 ) 



The problem of minimizing the drag given by equation (l 6 ), subject to the auxiliary 
n ^£y en ^ e Oration ( 19 ); can be expressed in the usual manner of the c al c ulu s of vi 


2A 


NACA TN 3189 






10 


NACA TN 3189 


An integration by parts, with f(-|3R2) = 0, yields 



, = A = Constant 

(i - W(i - n) 2 - p^i 2 


(23) 


The integral equation (23) is the same equation that appears in the 
minimum- induced-drag problem in lifting- line airfoil theory. 


■The solution of the integral equation 


is 



F(y)dy 

I - y 


= Constant 


F(y) 


Cy 4- Cl 

l/O - y)(y - a) 


Therefore, the solution of equation (23) is 


f(n) 


= (Ctj + Ci + P] 

» T1 + 


PRl - T) 


(24) 


q + PR2 

The condition that f(-f3R2) = 0 requires that Ci = C0R2J therefore, 


f(n) = 0/(11 + PR2K1 + P R i - 4) 


(23) 


The source-distribution function for .minimum wave drag is part of an 
ellipse. The constant C is fixed by satisfying the auxiliary condl 
tion, equation (19). The final result for the source-diBtribution 
function for minimum 'wave drag is 


f(n) = 


4u(ri 2 - r 2 2 ) 


rt(l - pRq + pR2)(l + 3£®i + PR 2 ) 


/(t + PF2)( 1 + P R 1 " 4) (26) 


Figure 2 which has an arbitrary vertical scale shows the source distribu- 
tion corresponding to the transition section of minimum wave drag. 



Minimum-Drag Shape 


The expression for the minimum-drag shape Is analogous to the auxiliary condition given by 
equation (17) except that a control surface of arbitrary radius R is used, where R is the 
radius of the section at x = X (see fig. l); the shape is given by 

pUrt^R 2 - R2^) — 


P(R-R2) 


2btpRv r (x,R)dx 


= 2j(p 


r 

nx-pR 

f'(i)(x - t)ae 

y p(R-R 2 ) 


/(x - if - P £r2 


(27) 


Interchanging the order of integration, integrating over x, and substituting the source- 
distribution function yield 


P 2 (r 2 - R 2 2 ) = 


8p 2 (R x 2 - R 2 2) 


«(1 - pRi + fffi 2 ) (l + 3 PRi + PR 2) -PR2 


/ 


«-pR i jgg ■ g)^ x _ |) g _ p %E 

]/U + (5R 2 )(l + f®! - I) 


d£ 


(28) 

Since all the radii in equation (28) are multiplied by (3, one calculation gives a family of 
shapes for which, at a given value of X, |3 R(X) = Constant. Thus the minimum-drag shape obeys 
the Gothert similarity rule (ref. 11 ). 


At the general point (X,R), the integral in equation (28) involves elliptic integrals and 
cannot be evaluated in elementary form. However, at one particular value of X, corresponding 
to the center of the source-distribution ellipse, the integral reduces to elementary functions. 
Thus, for 


x = 1 + PRl - 


p 
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^2 = (p^i 2 + p%*g 2 )(i + gRi + ffig) 2 - Qp^i 2 ^ 2 

2 |7l + pRi + PR 2 ) 2 - S^ 2 *! 2 + P 2 *^ 


( 29 ) 


ro 


Equation (28) may be solved by numerical- graphical methods. By a change of variable of 
the form 


t = 1 + ffil - P^2 + 1 + ffil + P R 2 cog 6 

5 2 2 

equation (28) is converted into the following form suitable for numerical integration: 


P 2(r2 . r 2 2) „ *<*1* - R g g )( 1 + I”1 + P R 2-) 
*(l - pRl + PR 2 )(1 + 51% + f®2) 


. 1 f i(X-pB)-(l+pRi-pR g ) 
[_ I+PR1+PR2 


B /(x - Lt-gl - --g g - 1 + Pi . f jgg eos 9 ) g . (30) 


The slope of the transition sectidn at the general point (X,R) is found by differentiating 
equation (28). The result is 


a(P) pc 

IT'® 


nX-pR 

(k±m±J-&2. s)(x - i)cL| 


pX-p 

t )<LS 1 

li-pR 2 1 

/(l + P 2 )(i + Pi - s) \/(x - l) 2 - p% 2 

u 

u 

-p 2 ] 

1 1 + Pg)(l + Pi - 6) l/(x - 6) 2 - p2R 2 


( 51 ) 
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For the general point, equation (3l) is not integrable in elementary 
form. However, at the downstream end of the section (x = 1, R = Rq) 
it may he integrated directly. The result is 



R=Rq 


/®q(l + £Rq - pR 2 )(l + f&l + 33^) 

/s^p^Rq 2 - fi^B.2^) + 3 R 1 + {JR 2 



(32) 


At the upstream end of the section, equation (31) may he evaluated hy 
a suitable limiting process. For example, if R = R 2 and X = e 
(where e is a small quantity), the integrals in equation (3l) may he 
evaluated if e is neglected in comparison with quantities of the 
order of unity. After the evaluation, e is allowed to vanish, and 
the result is 


(MMl 

La* J x=0 

r=r 2 


1 + 


^3^2(1 “ 3% + 3R2)(i + 33Rq + 3^) 
1/2/1 + £Rq + pR 2 (p2Ri 2 - 3 2 R 2 2 / 


1 


(33) 


For R 2 = 0 (the 
slope of the Mach 


R 1 f r 2- 


projectile tip), the front slope is equal to 
cone. The hack slope is different from zero 


the 

for 


Minimum Wave Drag 

The expression for the minimum wave drag is obtained hy inserting 
the source-distribution function (eq. ( 2 6 )) into the drag integral 
(eq. (16)). The expression for the drag of the transition section is 


&tpU 2 (Rq 2 - R 2 2 ) 2 

D = (3b) 

(l - 3 r i + 3 R 2)(l + 33^1 + 3 r 2^ 


and the drag coefficient (based on the area of the section projected on 
a plane normal to its axis) is , 

4 


= 


Rq 2 “ r 2 2 


(l - pRq + p R 2 )(l + 33Rl + 3 r 2^ 


( 35 ) 
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Reversibility 

The results for the reverse-flow case may be obtained formally 
from the preceding direct-flow case by interchanging Rq and R2. If, 
in addition, X is replaced by 1 - X, the X-station will be measured 
from the same end of the configuration. The transit ion- sect ion- shape 
expression (eq. (28)} can be shown to be an invariant relation between 
X and R under the transformation 


R l— *R2 
R 2 -»Ri 
X— »1 - X 

Since the present analysis is a linear treatment, the reversibility of 
the minimum-drag shape cannot be considered as established to an order 
higher than the first. 

The drag coefficient (eq. (35)) is not symmetrical in Rq and R2- 
However, p^Cp is symmetrical to the first order in the small quanti- 
ties pRq and pR 2 or, if pRq and PR2 are not small, p^D is 
symmetrical, to the first order in the quantities pRq - PR2 and 
pR 2 - pRq. The minimum wave drag thus is the same to the first order 
in the appropriate small quantities for the direct-flow and reverse- 
flow cases. 


DISCUSSION OF RESULTS 


The drag integral (eq. (12)) presented in this paper depends upon 
the source-distribution function in the same way (except for slightly 
different weights of the kernel functions) that the slender-body drag 
expression (eq. (15)) depends upon the cross-sectional-area derivative. 
Since bodies with area-derivative discontinuities can be generated by 
continuous source distributions, equation (12) gives finite drag for 
certain types of bodies (for example, bodies with shoulders forward of 
the base) for which the usual slender-body drag expression diverges. 

The drag integral presented herein is thus applicable to a larger class 
of bodies than the Ward drag expression. 

Some indication of the range of validity of equation (12) is given 
by a comparison of the values of drag calculated by various methods for 
several cones. Figures 3(&)> 5(b), and 3(c) present cone drag coeffi- 
cients plotted against pRq for semivertex cones of 5°> 10° , and 15°, 
respectively. The drags are calculated by the following methods: 
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Method A - The drag integral (eq. (12) ) vith mass- flow continuity 
(eq. ( 17 ) with R 2 = 0) used as the boundary condition. Since the con- 
tinuity condition is used in deriving the drag integral, this method 
uses a consistent boundary condition. 

Method B - The solution of the linearized supersonic-flow equa- 
tion for a cone (f ' ( |) = Constant) used with the exact boundary condi- 
tion { — = — — — I and the exact pressure equation. This method may be 

\dx U + v x / 

called the exact linear method. Van Dyke (ref. 12) has presented cone 
drags by this method. 

Method C - The exact characteristics method for conical flow as 
given by reference 13 . 

Method D - The drag integral (eq. (12) ) with the boundary condi- 
tion — = — — — . This method was used only for the 5° cone, 
dx U + v x 

Each of the figures also shows the drag of the minimum-drag projectile 
tip of the same thickness ratio as was determined in the present analysis 
(curve I in fig. 3 )« 

Figures 3(a), 3(b), and 3(c) show that the drag integral with the 
continuity boundary condition (method A) gives cone drags closer to 
the exact value (method C) than the exact linear method (method b) for 
values of gRq less than approximately 0.30. From a comparison of 
methods A and D for the 5 ° cone, the divergence of the c one drag calcu- 
lated by method A as pR^ — >1.0 is seen to be caused by the use of the 
mass-flow continuity boundary condition. These results indicate that 
equation ( 12 ) with equation ( 17 ) as the boundary condition gives good 
first-order results for Mhch numbers above the transonic range and less 
than the Mach number at which the tangent of the Mach angle is less 
than approximately three times the slope of the body, that is, a moderate 
supersonic range. 

Projectile tip (Rg = 0).- Figure k presents two families of minimum- 

drag projectile- tip shapes (pRq = 0.2 and O.t) calculated by equa- 
tion (30). The Von Karman projectile tip and the limiting case PRi = 1.0 
are shown for comparison. The Mach number and thickness-ratio dependence 
are related by the fact that the minimum-drag shape obeys the Gothert 
similarity rule (ref. ll). The slope at the base varies from zero to 
the cone slope as j3Rj_ varies from 0 to 1.0. 

In both the Von Karman analysis and the present analysis the basic 
assumptions of linearized theory are strongly violated at the nose of 
the projectile tip. Specific results of either analysis at the nose 
therefore cannot be accepted with confidence. However, inasmuch as 
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linear theory often qualitatively predicts correct trends even outside 
its range of quantitative validity, it is interesting to note that in 
both analyses the nose of the minimum- drag projectile tip is as blunt 
as the analysis will allow (infinite slope for the Von Karman tip and 
tangent to the Mach cone in the present analysis). It is probably 
true that the exact shape of a projectile tip for mini m u m wave drag 
has a blunt nose. 


Figure 5 presents the ratio of the m in imum- wave-drag coefficient 
based on the base area to the Von Karman drag coefficient as a function 
of pR]_. The ratio has a minimum value at pRq = 0.33 and diverges 

as fJRi — >1.0. The results for the cone discussed previously suggest 
that the drag diverges as pRq — >1.0 because of the use of mass-flow 
continuity (eq. (17)) for the boundary condition and that the results 
are reasonably good for pRq less than about 0.3- 


General case .- Figure 6 presents the radius at an intermediate 
point and the slopes at the two ends of the minimum- external- wave-drag 
transition section for four configurations with pRq - PR2 = 0.4. Fig- 
ure 7 presents these three quantities for five configurations for which 
pRq - PR2 = 0.05. Equations (29), (32), and (33) indicate that for 


pRl - PR2 < — and PRq + PR2 > 


(pRl - pR 2 )(l + pRi - pRg) 
1 - 3(fBi - PR2) 


the shape 


is a reverse curve lying within the truncated cone connecting the two 
cylinders for a distance downstream of the nose which increases as 
pRl - pRg decreases. The entire shape was calculated by equation (30) 
for the case pRq - PR2 =0.05 and pRq + PR2 = 0.118, which corre- 
sponds to the minimum front slope for pRl - PR2 = 0.05 (shown in 
fig. 7). This reverse-curve effect is more pronounced in the range 
of small values of pRq - PR2, the range where the linearized treatment 
is certainly correct. 


For a comparison with the general case, the drag of a truncated- 
cone section for pRq - pR2 = 0.05; P^l + P^2 = 0.1l8, and P = 1 was 
calculated by the method of reference 14. The truncated- cone drag 
coefficient (0.0228) is l8 percent higher than the minimum drag 
coefficient (0 . 0193 ) • 


CONCLUDING REMARKS 


The drag Integral for bodies of revolution at supersonic speeds is 
presented In a double-integral form depending explicitly upon the 
source-distribution function. The drag integral Is applicable to a 
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larger class of bodies than the usual slender-body drag expression and 
reduces to the Ward drag expression in the proper limit. Results for 
cones indicate that the drag integral with mass-flow continuity as the 
boundary condition gives good first-order results in the moderate 
supersonic range. 

The minimum- wave -drag problem for a transition section connecting 
two semi-infinite cylinders is solved with the aid of a suitably chosen 
auxiliary condition. The source distribution, the minim um drag, the 
slopes at the two ends of the section, and the radius at an intermediate 
point are obtained in terms of elementary functions . The entire shape 
is obtained in an integral form amenable to numerical evaluation. The 
minimum- drag shape obeys the Gothert similarity rule. The minimum-drag 
shape and the minimum drag are unchanged to the first order when the 
flow direction is reversed. The projectile tip is a special case and 
is compared with the Von Karman projectile tip. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 1, 195^* 
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APPENDIX 

INTEGRATION IN EQJJATION (ll) 


The Integral of equation (ll) is 


I = 


l+ £?l 1/^ 

Tl+!3Rl v 


(x - | + x - i})cbc 

|) 2 - pSRi 2 vtx - "-"^) 2 - p^q 2 


(Al) 


where the lower limit is the larger of the two quantities | + pRq or 
T) + pRq. The equation may he evaluated in the following manner. In 
order to convert the radicand in the denominator of the integrand into 
a quadratic function of the square of the integration variable, the 
variable is changed by using 


x = S - 


I + T] 
2 


Substitution of equation (A2) into equation (Al) gives 


(A2) 



where Sjy and Sq are the appropriate limits for the variable S. 
Since the numerator of the integrand in equation (A3) is an odd 
polynomial in S, considerable simplification is effected by taking a 
new variable z proportional to S 2 , such that 
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Substituting equation (Pk) into equation 



Ti+pRl 


(x - |)(x - q) - 
- Tl) 


1+pRl 
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Figure 2.- Source distribution (arbitrary vertical scale) corresponding 
to transition section of minimum drag. Ellipse of which the distribu- 
tion is a part is shown by dashed curve. 
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Drag coefficient 



(a) 5° cone. 


Figure 3-- Comparison of drag coefficients for three cones determined 
by several methods and for minimum- drag projectile tip of the same 
thickness ratio as the cone. Method A - Drag integral vith mass- 
flow- continuity boundary condition. Method B - Exact linear method. 
Method C - Exact characteristics method from reference 13. Method D - 

dr v r 

Drag integral with = jj- + . I - Drag coefficient of minimum- 

drag projectile tip. 
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Figure 7.- Radius at intermediate point and slopes at two ends of minimum- 

external-drag transition section for five configurations. pRi - pR 2 = 0.05. 
Entire contour shown for pRj_ + PR2 = 0.118. 
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